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We introduce a simple and predictive model for determining the phase stability of ternary 
phospholipid-cholesterol mixtures. Assuming that competition between the liquid and gel order of 
the phospholipids is the main driving force behind lipid segregation, we derive a Gibbs free-energy of 
mixing, based on the thermodynamic properties of the lipids main transition. A numerical approach 
wax devised that enable the fast and efficient determination of the ternary diagrams associated with 
our Gibbs free-energy. The computed phase coexistence diagram of DOPC/DPPC/cholesterol re- 
produces well known features for this system at 10° C, as well as its evolution with temperature. 
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Lipid membranes have captured the attention of physi- 
cists, chemists and biologists alike owing to their promi- 
nent structural and functional role in living cells, and to 
the pervading use of lipid based products in the pharma- 
ceutical, food and cosmetic industries In the recent 
years, ternary lipid systems comprising phospholipids 
and cholesterol were established as important model sys- 
tems in relation with lateral lipid segregation phenomena 
in membranes. At the root of this rich phenomenology 
lies a complex phase behavior of certain lipid mixtures, in 
which cholesterol molecules play a unique role. Ternary 
systems comprising a phospholipid in the liquid state, a 
second phospholipid in the gel state and cholesterol dis- 
play several unique characteristics, the two most impor- 
tant of them being arguably a wide binary coexistence 
gap between the so-called liquid disordered (Ld) and liq- 
uid ordered (Lq) phases, and a region of triple coexis- 
tence [3 5]. 

A prominent feature of phospholipid membranes is 
their main thermodynamic transition, from a low temper- 
ature gel or rippled gel phase (usually termed Lp^Pp), to 
a high temperature liquid phase (La) [I]- A sharp change 
in the statistical conformations of the hydrophobic chains 
occurs at the transition, associated with membrane thin- 
ning and an increase of the area per phospholipid. The 
transition takes place at a melting temperature Tm whose 
precise value depends on the chemical structure of the 
chains. It is a weakly first order transition for mem- 
brane of pure composition. The thermotropic behavior 
of phospholipid binary mixtures have been studied since 
the seminal work of Mabrey and Sturtevant @ and sup- 
ports the view that liquid-prone and gel-prone lipids tend 
to phase separate. This suggests that mixtures composed 
of phospholipids with same head group have their mis- 
cibility behavior governed by the competition between 
their liquid and gel conformations. Hence, a description 
of ternary lipid systems must necessarily, and to some ex- 
tent, account for the internal state of the phospholipids. 

We introduce in this Letter a simple, robust and pre- 



dictive model of phospholipid-cholesterol mixtures phase 
coexistence. Our model aims at predicting the misci- 
bility behavior of phospholipids, based upon the ther- 
modynamic properties of their main transition 
The thermotropic behavior of the phospholipid tails is 
modeled with an internal order parameter, in a way 
that schematically reproduces the thermodynamic fea- 
tures of the liquid-gel transition. The peculiar features 
of cholesterol-phospholipid interactions [l2| are ac- 
counted for by means of a non-ideal mixing entr opy gen- 
eralizing the one proposed by Komura et al. [l3|. In 
addition, we outline an original numerical procedure for 
the derivation of ternary phase diagrams. 

The stability of a mixture is governed by its Gibbs free- 
energy of mixing Gmix- This thermodynamic function 
expresses the free-energy difference between the homo- 
geneous mixed state and the system comprising the sep- 
arated pure components with same global composition. 
Ternary mixtures in membranes are therefore described 
by Ginix(T, ni,n2, n^), where the Ui are the molar content 
of each lipid, and T the temperature. We assume that the 
hydrostatic pressure remains constant, close to the stan- 
dard value, and that lipid membranes have no tension, 
so that these two intensive parameters are absent from 
Gmix- Membranes are assumed to be well hydrated with 
a large excess of water, and free from any other external 
fields ^14] . 

The simplest realization of our model, described in this 
Letter, assumes that the areas per lipid do not change 
with temperature. This allows us to express Gmix in 
terms of the surface fractions , (f)2 and ^c" of each lipid 
component, 4>c' standing for the cholesterol. We claim 
that it is acceptable to neglect all changes in the area 
per lipid at T = T^, in a first approximation to the 
description of the mixture phase coexistence, especially 
in the absence of membrane tension. This restriction is 
not a fundamental one, and versions of our model ac- 
counting for membrane expansion at the transition will 
be presented elsewhere [3]. 
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A frequent approximation to non-ideal solutions is ob- 
tained by combining an ideal mixing entropy with a 
quadratic enthalpic term. For a membrane of area A, the 
free-energy of mixing per unit of area, g^i^ — /SGnii^/A, 
reads: 



^ ln(0i) + ^ ln(02) + — ln(0c) 
ai a2 ac 

+Xi2 4'i4>2 + Xic Mc + X2C 4>24>c , (1) 



with 4>i+4'2+(t'C = 1 and 13 the Boltzmann factor l/ksT, 
and ai , 02 and ac are the areas of the three different lipid 
species. One observes that expression ([!]) alone does not 
account for the main transition, and that the effective 
interaction parameters Xij are expected to display a non 
trivial and unknown temperature dependence. We could 
not obtain, based on ([T]), any diagram comparable with 
the experimental ones. Moreover, our attempts to do 
so lead to unexpectedly large values of ac, conflicting 
with the view that cholesterol molecules are smaller than 
phospholipids. 

Experiments and simulations either suggest that 
cholesterol may form transient bound complexes with 
phospholipid moieties surrounding them [ill [l2j . The 
presence of such complexes prevents cholesterol molecules 
from occupying too many neighboring positions, thus giv- 
ing the semblance of a large excluded volume. To ac- 
count for the statistical behavior of cholesterol molecules, 
we derived the following entropic contribution, where all 
lipid areas have been set to unit: 



G^""^ [4^1 A2) = 4^1 H4i) + h H42) + 0C7 ln(0c) 
- (^1 + 02) l0g((/)l + (j)2) + (<^"^ - 0c) log(l - 5cj)c). 



(2) 



This expression originates from a lattice enumeration of 
configurations, following Flory, where each cholesterol 
molecule deprives an average number 5 of neighboring 
sites from the possibility of being occupied by another 
cholesterol molecule. Expression ([2]) reduces to ideal mix- 
ing when (5^-1, and is similar to the entropy of Komura 
et al. for binary systems and 5 = 2 [l3| . It accounts 
for both the mixing term and the excluded volume be- 
tween cholesterol molecules, consistent with the presence 
of cholesterol complexes. It finally sets a maximal value 
<f>c < 1/5 to the mixture content in cholesterol. 

Several Landau order parameters have been proposed 
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for describing the liquid-gel chain melting [7|, y 
[lij . Theories based on a single scalar order parameter 
can pretend to capture only rough features of the main 
transition. This is why, for simplicity and generality, we 
restrict ourselves to a schematic approach, a two-state 
model based on a scalar order parameter m restricted to 
the interval [—1,1]. The mean-field Ising model under 
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FIG. 1: (Color online) (a) Discrete approximation to the 
ternary phase diagram obtained from the lower convex hull 
method explained in the text, (b) Final phase diagram deter- 
mined from the analytical refinement of the discrete approxi- 
mation. 



finite external field serves as a guide [17|. 
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hi{T) 



All AM 

-^{T~T,);h2{T) = ^{T-T2), (4) 



and R is the gas constant. Eq. ([3]) represents a mean- 
field Ising model in the presence of an effective external 
field hes = hi(t>i + /i2(/'2- For J > J^. = 1/4, the m 
dependence of Q^^^ is shaped as a double well potential, 
with a linear bias induced by hcS- The order parameter 
m fluctuates freely around an optimal value m* located 
at the global minimum of G^^\ with (T, 0i, 02) constant. 
In practice, as soon as J is larger than 0.3, the two lo- 
cal minima are close to ±1. We conventionally ascribe 
rUg ~ —1 to represent the gel phase, while to/ ~ 1 is 
set to represent the La liquid phase. The coupling J ex- 
presses that neighboring lipid tails tend to be in the same 
state, irrespective of the chemical species to which they 
belong. This is the main cooperative effect that induces 
the liquid-gel separation. 

The term /igff acts as a switch that drives the order pa- 
rameter from its m* = m,g gel state to its to* = to/ liquid 
state, upon changing its sign from negative to positive. 
One sees from @ that if lipid 1 is pure, hi (T) changes its 
sign at T = Ti while m*{T) suddenly jumps from m,g to 
TO/, with an amplitude Ato* ~ 2. The enthalpic change 
is close to A77i, so that A7?i can be identified with the 
molar latent heat at the main transition. The difference 
in Gibbs free-energy between the liquid and gel phases 
at T 7^ Ti is then extrapolated linearly with T at the 
vicinity of the transition. The contribution C/'^' is thus 
entirely specified by the temperature and the enthalpy 
difference at the transition of the lipid species composing 
the mixture. This term induces a liquid-gel separation in 
the whole temperature range [Ti , ■ 
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Cholesterol-phospholipid interactions are further ac- 
counted for with a specific interaction term 

C/(<^)(0i,02,™) = xMc - + (l>2)4>c (5) 

Eq. ^ comprises a m-dependent contribution, propor- 
tional to ^, which marks the preference of cholesterol 
for either the m — 1 liquid order, or the m = — 1 
gel order, according to the sign of ^. This interaction 
does not discriminate between saturated and unsatu- 
rated lipids, since it describes a coupling to the chain 
order. Other cholesterol-phospholipid interactions are 
accounted for by a m-independent contribution in x- 
This order-insensitive interaction introduces a coupling 
between cholesterol and the unsaturated component. We 
found that interactions between cholesterol and the satu- 
rated component are not required to describe the exper- 
imental phase diagrams. 

The sum ^(T, 01, 02, m) = + gt^) + G'^^^ is our 
main thermodynamic potential. In the absence of sur- 
face tension, and assuming that changes in internal state 
occur at constant volume, m fluctuates around its opti- 
mal value m*, with m ~ m* in the thermodynamic limit. 
Therefore, the free-energy of mixing (/mix (7^, 0i, ^2) is ob- 
tained by minimizing Q (T, 0i , 02 , m) with respect to m 
at fixed temperature and surface fractions [loj . 

gmix(r,0i,02) = niin [C/(T,0i,02,to)]It,0i,02 

= g(r,0i,02,™*(T,0i,02)). (6) 

The Gibbs free-energy surface ^mix arises from project- 
ing out m from the four-dimensional potential Q. Equiv- 
alently, ^mix can be viewed as the restriction of Q to 
the submanifold dmG = 0. For J > Jc, the constraint 
dmG = defines two stable submanifolds, a negative 
mg{T, 01, 02) and a positive mi{T, 0i, 02) branch, linked 
respectively with the gel and liquid phases. At fixed 
(T, 01,02), only one of the two branches is the stable 
one, while the other branch remains metastable. 

In ternary mixtures it is known that phase coexistence 
is dictated by the convexity properties of gmixiT, 0i, 02) 
relative to 0i , 02 . Minimization of ^mix is equivalent to 
finding ^mix, the lower convex hull of f/mix, that displays 
regions of three kinds. Points of the convex hull in con- 
tact with the original surface are stable monophasic re- 
gions. Regions lying above triple tangent planes corre- 
sponds to triple coexistence. The remaining part of the 
convex hull correspond to developable patches of surface, 
wrapped around double tangent lines that connects pairs 
of points with distinct composition (tie- lines) In 
practice, we used a public domain routine, qhull, to com- 
pute a discrete approximation of ^mix [l9[ . The starting 
point is a fine mesh discretization of gmix over the relevant 
domain of composition {0i, 02 > 0, 0i -I- 02 < 1 — 1/5}. 
Then qhull computes a triangulated surface approach- 
ing ^mix- Triple coexistence correspond to facets with 
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FIG. 2: (Color online) Computed evolution of the 
DPPC/DOPC/Chol triangular diagram with temperature. 
The computed diagrams display the main expected features 
for this lipid mixtures. 



all sides much larger than the initial mesh size. Double 
coexistence is associated with elongated triangles, with a 
shortest side much smaller than the two longer sides, the 
latter being oriented parallel with the tie-lines. Finally, 
small facets of the convex minimization are linked with 
stable monophasic regions. The projection of the trian- 
gulated gmix surface onto the composition plane provides 
directly a fast and accurate approximation of the Gibbs 
phase diagram shown in fig. [l](a). This discrete solu- 
tion can then be refined according to the usual rules for 
the determination of ternary phase diagrams. Details of 
these calculations will be given in a forthcoming publica- 
tion 

Arguably, the most widely studied ternary lipid phase 
diagram for membranes containing cholesterol is the 
ternary mixture DOPC/DPPC/Chol. DOPC is a dou- 
ble tail molecule with 18 carbons per tail and a cis- 
unsaturated bond on the ninth carbon. It has a low main 
transition temperature Ti at — 21°C and an associated 
enthalpy change AHi evaluated at 7.7 kcal/mol (20| . 
DPPC has also a phosphatidylcholine head but here the 
two 16-carbon tails have no unsaturated bonds, leading 
to a higher transition temperature T2 at 42° C, and asso- 
ciated enthalpy change AH2 = 8.7 kcal/mol I20'|. Within 
our approach, the above temperature and enthalpy val- 
ues and a single coupling parameter J — 0.35 - see eq. ([3]) 
- describe well published data 2l| for the experimental 
binary DOPC/DPPC mixing behavior. For a complete 
description of the ternary phase diagram our model fur- 
ther requires three parameters: (5,^ and The value of 
6 is bound by the maximum amount of cholesterol that 
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a membrane can hold, estimated here at 4>c — 0.45 or 
equivalently (5 = l/0c = 2.2 j|. The parameter ^, as 
defined in eq.jS]), quantifies the preference of cholesterol 
for the liquid disordered phase with respect to the gel 
phase [l|,[2§|. Binary diagrams of DPPC/Chol computed 
with ^ = 1.8 reproduce well the temperature behavior of 
the miscibility gaps occurring between the i/j-gel phase 
and the Lq liquid ordered phase as determined from Gi- 
ant Unilamellar Vesicles in [23| and from NMR 2J| - see 
fig. ([2]). Eq.© also includes a repulsive m-independent 
interaction between cholesterol and the unsaturated com- 
ponent, given by x- At 10°C a value of x = 6.0 pro- 
motes both a triple coexistence region and a large bipha- 
sic domain, with oblique LojLd tie lines of positive slope 
and a plait point on the left side of the Gibbs triangle, 
in agreement with experimental findings 's']. A larger 
value of X would open there a miscibility gap, unseen 
for DOPC/Chol, while a smaller value would drastically 
reduce the biphasic Lq/ Ld domain. In our model a con- 
tinuous path connects the Lo and Ld regions around the 
critical point, where the order parameter m remains posi- 
tive. The Lfj phase on the lower right corner of the phase 
diagram is associated with a negative m-value and sepa- 
rated from the rest of the diagram by a discontinuity in 
m across the coexistence gap. 

The evolution of the Gibbs diagram with respect to 
temperature is displayed in fig. The diagram was 

drawn here under the assumption of a constant ^ param- 
eter, and a constant value of the coupling contribution 
ksT ■ X denoting its pure enthalpic nature. Our diagram 
at the lower temperature in fig. ([2]), displays the main 
expected features for a DOPC/DPPC/Chol mixture. As 
the temperature increases, the triple coexistence trian- 
gle shrinks and vanishes between 15°C and 20°C. Also, 
the Lo/Ld biphasic region, resting initially on the left 
side of the triangle, detaches and assumes a closed shape 
of decreasing surface. A first critical point, on the left 
side of the Gibbs diagram, evolves slowly away from the 
edge (3]. A second critical point, initially hidden by the 
triple coexistence region, is predicted to emerge at the bi- 
nary detachment temperature. Note also that within our 
approach the extension of the Lp domain, on the right 
bottom part of the diagram, gradually decreases up to 

As a summary we presented a model for phase coexis- 
tence in ternary phospholipid-cholesterol mixtures. Our 
approach is based on a mean-field two-state description 
of the phospholipids and a non-ideal mixing entropy for 
cholesterol. All but one parameter in our system are 
bound by known thermodynamic or physical properties. 
We combined analytical theory and numerical convex 



minimization to show the predictive power of our model 
for the ternary mixture DOPC/DPPC/Chol. Resulting 
phase coexistence diagrams reproduce well known fea- 
tures for this system at 10°C as well as the temperature 
variation of the diagram. We are also encouraged by pre- 
liminary results from related ternary systems that further 
indicate that our model successfully describes a large va- 
riety of lipid/lipid/cholesterol mixtures. 

*fabrice. thalmann@ics-cnrs.unistra.fr 
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